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56 CO 7-RAYS FROM SN1998BU: COMPTEL UPPER LIMITS 
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ABSTRACT 



The type la supernova SN1998bu in M96 was ob- 
served by COMPTEL for a total of 88 days starting 



17 Hays after the first detection of the supernova. 



The accumulated effective observation time was 4.14 
Msec. The COMPTEL observations were performed 
in a special instrument mode improving the low- 
energy sensitivity. We generated images in the 847 
kcV and 1238 keV lines of 56 Co, using improved point 



sproad functions for the low-energy mode. We do not 
detect SlNiyy8bu. A spectral analysis ot our data 
also confirms the non-detection of the supernova. We 
discuss the event for which our upper limits constrain 
the standard supernova models. 
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1. THE OBSERVATIONS 



On May 9.9 UT 1998 (TJD 10942.9) [M. Villi (1998| ) 



obs erved a supernova in M96 (NGC 3368), which 
wa s labelled SN1998bu. — Fiom sp cctiogianis it w ss 



classi fied as type la supernova ( Ayani fc Meikle 
1998). Based on pre-explosion observations and an 



estimate of the m aximum blue light at tn 



10952.7 ± 0.5TJD Meikle & Hernandez (1 



de- 
termined the explosion dat e to be May 2 ± 1 UT 



1998 (i.e. TJD 10935 ±1). [Hjorth fc Tanvir (1997Q 
determined the distance to M96 by Cephcid mea- 
surements to 11.3 ± 0.9 Mpc. 

Observations of supernovae in the optical and neigh- 
bouring wavelength bands concentrate on informa- 
tion on the light curves from such events. The time 
evolution of theses light curves could be understood 
as being powered by reprocessed 7-ray line emission 
from freshly synthesised short-lived radio-isotopes 
(e.g. decay chain of 56 Ni). However, due to the 
creation of the low energy photons by secondary pro- 
cesses long after the explosion itself most of the infor- 
mation on the initial state right after the explosion is 



lost. Therefore the distinction of different supernova 
models via their predicted light curves alone is uncer- 
tain. In contrast 7-ray lines carry information from 
the very early phase of the supernova and m ay allow 
a discrimina tion of the theoretical models (Hoflich 
1998). Depending on the particular type la 



et al. 



supernova model and the sensitivity of current days 
7-ray instruments (OSSE & COMPTEL) the obser- 
vations of 7-ray lines are limited to a maximum dis- 
tance of ~ 15 Mpc. Up to now only one type la su- 
pernova (SN19 91T) was marginally detected (Mor- 
ris et al., 1997|). According to the distance estimate 



SN1998bu opened a second chance to observe low en- 
ergy 7-ray lines from a type la supernova. 
The COMPTEL observations in direction of M96 
started 17 days after the explosion (TJD 10952) and 
covered the time until TJD 11071 (i.e. 136 days 
after the SN). Due to some breaks we had in to- 
tal 88 days of supernova observations, which sums 
to an effective observation time of ~ 4.14 • 10 6 s. 
To increase COMPTELs low energy sensitivity the 
telescope was switched into an dedicated low-energy 
mode, decreasing the module thresholds of the D2 
detectors well below 600 keV. (For a deta iled descrip- 
tion of the COMPTEL instrument see [Bchonfcldor 
et al. (1993p ). Due to the late start of the observa- 
tion program COMPTEL missed the decay of 56 Ni 
(77/2 = 8.8 d). So the observations were focused on 
the detection of 847 a nd 1238 keV lines from the sec- 
ondary isotope 56 Co ( Gomcz-Gomar et al., 19981 ) . 



2. IMAGING ANALYSIS 



The observations have been performed in a special 
'low energy' mode, in which the D2 module thresh- 
olds were lowered to the minimum. By using ap- 
propriate point spread functions taking care of the 
real hardware thresholds we gain in low energy sen- 
sitivity. Due to the strong contamination of the D2 
data due to the 511 keV background line in some D2 
modules a higher low energy cutoff (560 to 600 keV) 
was applied for these modules. However, in imaging 
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Figure 2. Maximum-Likelihood images (left: 847 keV /right: 1238 keV) of SN1998bu source region in a local 
coordinate system centered on the direction towards SN1998bu. The SN position is marked with a star. The 
contour levels give the flux in units of 10 -8 phcm~ 2 s _1 . 2900 in the left map correspond to la, whereas 1600 
do so in the right. 



analysis another possibility is given by an adequate 
Compton scatter angle selection. Figure p] shows 
the different possibilities of 511 keV background 
suppression by cutting the data-space in a E1-E2- 
representation (For a description of the meaurement 
principle and its realisation see Schonfelder et al 
(19^3^). The band of contour-lines represents the 
distribution of observed events with a total energy 
deposit (sum of E-y and E 2 ) of 847 ± 75.3 keV. 
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Figure 1. Cuts in the energy data-space in the case 
of the 847 keV line. 



The 511 keV contamination could be clearly seen as 
an area of increased density below the horizontal 600 
keV line. This background could be suppressed ei- 
ther by an energy selection in E% or an cut in the 



scatter angle distribution tp. To really make use of 
the low mode data new point spread functions (PSF) 
have to be generated. The COMPTEL team to- 
day uses two approaches to generate point spread 
functions, one based on a Monte Carlo simulation 
of the detector response (SIMPSFs) and the other 
based on model assumptions using measured data 
(MODPSFs). In standard case with sufficient statis- 
tics both methods give rather identical point spread 
functions. However, we used of both methods to gen- 
erate low mode PSFs. 

We tested the adequacy of our response treatment 
on Crab observations taken in the same instrument 
mode. In the energy window of the low energy line 
we clearly detect Crab with both PSFs. However, 
due to a more accurate treatment of the detector 
settings in the response simulation, the image ob- 
tained with the SIMPSF is smoother than for the 
MODPSF. The point-source significance rises from 
3.02cr using a standard PSF to 5.1er by applying the 
new SIMPSF (for the MODPSF the value is 4.9cr). 
Also in the case of the high energy line the new PSFs 
increase the significance, with a smaller increase be- 
cause of the E~ 2 spectrum of Crab and the smaller 
effect of lower detector thresholds in the higher en- 
ergy regime, the significance rises from 0.63er to 1.6a. 
We applied both sets of point spread functions to the 
accumulated SN1998bu data using a maximum like- 
lihood reconstruction. Figure ^| shows the images 
obtained with the SIMPSFs. These images as well 
as the images obtained by applying the MODPSFs 
show no hint of a supernova detection. Also an addi- 
tional attempt using a combined point spread func- 
tion for the search of both lines simultaneously fails 
in detecting SN1998bu. We deduce 2a upper limits: 
From the flux distribution of all independent pixels 
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of t hese maps their variance is determined, assuming 



mate of M9 6 supports the deduced distance of Manoz 



a C aussian distribution. Using a Bayesian method 
( Gcorgii et al., 1997[ ), which accounts for systematic 



as well as statistical uncertainties, 2a upper limits 
of 3.7 • lO^phcm^s- 1 for the 847 keV line and 
1.9-lCr 5 phcm~ 2 s" 1 for the 1238 keV line are found. 



3. SPECTRAL ANALYSIS 



In addition to the imaging analysis we also performed 
a spectral analysis of the obtained data from the area 
around SN1998bu. In contrast to the imaging anal- 
ysis it is possible in the spectral domain to apply a 
(^-selection to suppress the 511 keV contamination. 
To achieve an adequate suppression we applied a en- 
ergy threshold to the D2 data at 600 keV. We then 
performed a spectral scan on a 3° wide grid around 
the position of SN1998bu. In each spectral analysis 
of this scan we used data from a 3° cone as source 
spectrum and data from a 3° — 7° cone mantle as 
background spectrum. After fitting the background 
spectrum to the source spectrum the residual spec- 
trum should contain any excess source signal. Subse- 
quently, we fitted template Gaussians, with a width 
corresponding to the instrumental energy resolution, 
to the residual spectrum to obtain the line intensities 
of the expected gamma-ray lines. Using this method 
no significant signals could be detected. Moreover, 
no significant differences between the source position, 
position of SN1998bu, and the off-source positions in 
scanning grid could be detected. 
We derived 2a upper limits by the means of the fol- 
lowing procedure: We generated histograms of the 
fitted intensities from all positions for both lines. 
The width of the distributions was then interpreted 
as a measure of the statistical and the systemati- 
cal uncertainties of the method. In that way we de- 
rived upper limits of 4.1 • 10 -5 phcm~ 2 s _1 for the 
847 keV line and 2.3 • 10~ 5 phcm" 2 s" 1 for the 1238 
keV line. The spectroscopically deduced upper limits 
are somewhat higher because of the lower sensitivity 
du e to the hard D2 energy cut. The (,5-cut in imag- 



ing | analysis increases the sensitivity a little because 



of the l/(-Ei + £ , 2)-dependence of the (^-distribution. 
In addition the background treatment in the imaging 
analysis is more reliable than in the spectral domain. 



4. INTERPRETATION AND CONCLUSION 



For a comparison of our upper limits with recent 
predictions of theoretical type la supernova models 
(Gomcz-Gomar et al., 1998; Kumagai, 1998) the dis- 
tance estimate to the SN is essential. Tanvir et al 

the host galaxy 



p995l ) 



estimated the distance to 
M96 by means of an HST Cepheid observation to 
11.6 ± 0.9 Mpc, which was revised by Hjorth fc Tan 
yir (1997|) to 11.3 ± 0.9 Mpc. However, a distance 
determination based on Planetary Nebulae observa- 
tion suggests a much clo ser distance of 9.6 ± 0.6 Mpc 



et al. (1999 ) based on a revised Cepheid distance cal- 



ibration. 



847 keV 




10950 



11000 



11150 



1238 keV 



6 

■r 5 

E 4 

O 
-C 

Q- 

•? 3 
o 

I 2 




W7 (Kumagai) — 




W7DT (Kumagai) - 




HeCD (Kumagai) - 




WDD2 (Kumagai) — 




Deflag. (Isern) o— 




Del. Det. (Isern) o- 


-0 


Det. (Isern) o- 




SubCH. (Isern) o- 


-o 


2o upper limit: 




Imaging — 




Spectral 





11050 11100 
TJD 



Figure 3. Predicted 847 and 1238 keV light curves 
from various explosion models. In addition the mod- 
els use two different radiation transfer prescriptions 
(Gomez-Gomar et al., 1998; Kumagai, 1998). The 
COMPTEL upper limits are plotted as horizontal 
lines, where the length gives the time coverage of the 
COMPTEL observations. 

We use the greater distance estimate of 11.3 Mpc for 
the comparison. Figure || shows the predicted light 
curves from different type la supernov a and radiation 
transfer mod els. As pointed out by Gomez-Gomar 



(Fcldmcicr et al., 1997). This closer distance esti- 



ct al. (1998) the different treatments of the radia- 
tion transfer problem through the expanding explo- 
sion shells is a major uncertainty in comparing the 
explosion models. Nevertheless we exclude all but 
the Sub-Chandrasekhar and deflagration models. In 
the case of a closer distance our results become even 
more constraining. 

In summary our measurement favours the deflagra- 
tion or Sub-Chandrasekhar models for SN1998bu 
and renders the HeCD model as rather improbable. 
For detonation models a smaller mixing than applied 
in tested models is needed to be compatible with our 
analysis. 
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